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Abstract 

We present a stability study of the class of multivariate self-excited Hawkes point processes, that 
can model natural and social systems, including earthquakes, epileptic seizures and the dynamics 
of neuron assemblies, bursts of exchanges in social communities, interactions between Internet 
bloggers, bank network fragility and cascading of failures, national sovereign default contagion, and 
so on. We present the general theory of multivariate generating functions to derive the number of 
events over all generations of various types that are triggered by a mother event of a given type. 
We obtain the stability domains of various systems, as a function of the topological structure of the 
mutual excitations across different event types. We find that mutual triggering tends to provide a 
significant extension of the stability (or subcritical) domain compared with the case where event 
types are decoupled, that is, when an event of a given type can only trigger events of the same 
type. 
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I. INTRODUCTION 



Many natural and social systems are punctuated by short-lived events that play a par- 
ticularly important role in their organization. Such events can be conveniently modeled 



mathematically by so-called point processes [l, 2|. They are also called shot noise in physics 
3 s|| or jump processes in finance and in economics jo]. These models are characterized by 
their (conditional) rate X(t\H t ) (also called "conditional intensity") defined as the limit for 
small time intervals A of the probability that an event occurs between t and t + A, given 
the whole past history H t . In mathematical notations, this reads 

X(t\H t ) = liniA^o — Pr(event occurs in [t,t + A]\H t ) , (1) 

where Pr(X\H t ) represents the probability that event X occurs, conditional on the past 
history H t . The symbol H t represents the entire history up to time t, which includes all 
previous events. This definition is straightforward to generalize for space-dependent inten- 
sities \(t,r\H t ) and to include marks such as amplitudes or magnitudes (see below). The 
standard Poisson memoryless process is the special case such that \(t\H t ) is constant, i.e., 
independent of the past history. Clustered point processes generalize the Poisson process 
by assuming that the series of events are generated from a cluster center process, which is 
often a renewal process, and a cluster member process. 

The class of point processes that we study here was introduced by Hawkes in 1971 {?]- 



10j . It is much richer and relevant to most natural and social systems, because it describes 



"self-excited" processes. This term means that the past events have the ability to trigger 
future events, i.e., X(t\H t ) is a function of past events, being therefore non-mar kovian. Many 
works have been performed to characterize the statistical and dynamical properties of this 
class of models, with applications ranging from geophysical 11- 16|. medical [17] to financial 
systems, with applications to Value-at-risk modeling [isj], high-frequency price processes 



191 ] . portfolio credit risks 



and yield curve dynamics 



23]. 



cascades of corporate defaults 



2l|, financial contagion 22], 



While surprisingly rich and powerful in explaining empirical observations in a variety of 
systems, most previous studies have used mono-variate self-excited point processes, i.e., they 
have assumed the existence of only a single type of events, all the events presenting some 
ability to trigger events of the same type. However, in reality, in many systems, events come 



2 



in different types with possibly different properties, while keeping a degree of mutual inter- 
excitations. Among others, this applies to geo-tectonic deformations and earthquakes, to 
neuronal excitations in the brain, to financial volatility bursts in different assets, to defaults 
on debts in some firms or some industrial sectors, to sovereign risks in some countries within 
a currency block, to the heterogeneity of activity of bloggers on the Internet, and so on. 

These observations suggest that multivariate self-excited point processes, which extend 
the class of mono-variate self-excited point processes, provide a very important class of 
models to describe the self-excitation (or intra-triggering) as well as the mutual influences 
or triggering between different types of events that occur in many natural and social systems. 
These considerations have motivated us to present recently the first exact analysis of some 
of the temporal properties of multivariate self-excited Hawkes conditional Poisson processes 
2^ |. as they constitute powerful representations of a large variety of systems with bursty 
events, for which past activity triggers future activity. The term "multivariate" refers here 
to the property that events come in different types, with possibly different intra and inter- 
triggering abilities. Ref. 24j was a first step towards a systematic study of the multivariate 



self-excited point processes, first mentioned by Haw! 



KCS 
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limself in his first paper |7j], whose 



25] 



full relevance has only been recently appreciated 

The present paper is a complementary study to our previous paper [24(, which was fo- 
cused on temporal properties, by studying the general stability conditions of this class of 
models. Section 2 recalls the definition and notation of Hawkes processes, starting from 
the monovariate version and extending to the general multivariate formulation. Section 3 
presents the formalism of multivariate generating functions to derive the number of events 
over all generations of various types that are triggered by a mother event of a given type. 
Section 4 gives the stability conditions using the mean numbers of events of all generations. 
Subsection 4A provides the general relations. Subsection 4B studies the case of symmetric 
mutual excitation abilities between events of different types. Subsection 4C restricts to the 
case of just two different types of events, that allows an in-depth analysis of the new features 
resulting from the inter-type excitations. Subsection 4D presents the results obtained for 
a one-dimensional chain of directed triggering in the space of event types. Subsection 4E 
generalizes subsection 4D by studying a one-dimensional chain in the space of event types 
with nearest-neighbor triggering. Subsection 4F presents a quantitative measure of the size 
of the subcritical domain that allows us to study the influence of the inter-type coupling 
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strength. Section 5 concludes. 



II. DEFINITIONS AND NOTATIONS FOR THE MULTIVARIATE HAWKES 
PROCESSES 

A. Monovariate Hawkes processes 

Self-excited conditional Poisson processes generalize the cluster models by allowing each 
event, including cluster members, i.e., aftershocks, to trigger their own events according to 
some memory kernel h(t — ti). 

x(t\H t ,e) = x c (t) + ^2h(t-u,) , (2) 

i\ti<t 

where the history H t = {U}i<i<i u t it <t<u t+1 includes all events that occurred before the 
present time t and the sum in expression tfZfj runs over all past triggered events. The set of 
parameters is denoted by the symbol 0. The term A c (i) means that there are some external 
background sources occurring according to a Poisson process with intensity A c (i), which may 
be a function of time, but all other events can be both triggered by previous events and can 
themselves trigger their offsprings. This gives rise to the existence of many generations of 
events. 

Introducing "marks" or characteristics for each event leads to a first multidimensional 
extension of the self-excited process ([2]). The generalization consists in associating with each 
event some marks (possible multiple traits), drawn from some distribution p(m), usually 
chosen invariant as a function of time: 

X(t, M\H t , 0) = p(M) I A c (t) + £>(*- U, , (3) 

\ i\U<t J 

where the mark Mj of a given previous event now controls the shape and properties of the 
triggering kernel describing the future offsprings of that event i. The history now consists 
in the set of occurrence times of each triggered event and their marks: H t = {t iy Mj}i<j<jv- 
The first factor p(M) in the r.h.s. of expression writes that the marks of triggered 
events are drawn from the distribution p(M), independently of their generation and waiting 
times. This is a simplifying specification, which can be relaxed. Inclusion of spatial kernel 
to describe how distance impacts triggering efficiency is straightforward. 
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From a theoretical point of view, the Hawkes models with marks has been studied in 
essentially two directions: (i) statistical estimations of its parameters with corresponding 



residual analysis as goodness of fits |39l448j ; (ii) statistical properties of its space-time dy- 



namics 
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The advantage of the self-excited conditional Hawkes process includes a very parsimonious 
description of the complex spatio-temporal organization of systems characterized by self- 
excitation of "bursty" events, without the need to invoke ingredients other than the generally 
well-documented stylized facts on the distribution of event sizes, the temporal "Omori law" 
for the waiting time before excitation of a new event and the productivity law controlling 
the number of triggered events per initiator. 

Self-excited models of point processes with additive structure of their intensity on past 
events [lo| make them part of the general family of branching processes 49]. The crucial 
parameter is then the branching ratio n, defined as the mean number of events of first 
generation triggered per event. Depending on applications, the branching ratio n can vary 
with time, from location to location and from type to type (as we shall see below for the 
multivariate generalization). The branching ratio provides a diagnostic of the susceptibility 
of the system to trigger activity in the presence of some exogenous nucleating events. 

We refer in particular to Ref. [l^] for a short review of the main results concerning the 
statistical properties of the space-time dynamics of self-excited marked Hawkes conditional 
Poisson processes. 



B. Multivariate Hawkes processes 

The Multivariate Hawkes Process generalizes expressions fl3]) into the following general 
form for the conditional Poisson intensity for an event of type j among a set of m possible 



types (see the document 26j for an extensive review): 

Xj(t\H t ) = A°(t) + / M* - s ) 9k(x) N k (ds x dx) , (4) 

k=1 J (-oo,t)x7J 

where H t denotes the whole past history, A° is the rate of spontaneous (exogenous) events of 
type j, sources of immigrants of type j, A k j is the (k, j)'s element of the matrix of coupling 
between the different types which quantifies the ability of a type fc-event to trigger a type j- 
event. Specifically, the value of an element Aj^ is just the average number of first-generation 



events of type j triggered by an event of type k. The memory kernel hj(t — s) gives the 
probability that an event of type k that occurred at time s < t will trigger an event of 
type j at time t. The function hj(t — s) is nothing but the distribution of waiting times 
(here between the impulse of event k which impacted the system at time s, the system 
taking a certain time t — s to react with an event of type j, this time being a random 
variable distributed according to the function hj(t — s). The fertility (or productivity) law 
gk{x) of events of type k with mark x quantifies the total average number of first-generation 
events of any type triggered by an event of type k. We have used the standard notation 

I(-oo,t)xnf^ x ) N ( ds x dx ) '■= Efc|t h<t /(*i^i)- 

The matrix A^- embodies both the topology of the network of interactions between dif- 
ferent types, and the coupling strength between elements. In particular, includes the 
information contained on the adjacency matrix of the underlying network. Analogous to the 
condition n < 1 (subcritical regime) for the stability and stationarity of the monovariate 
Hawkes process, the condition for the existence and stationarity of the process defined by 
01]) is that the spectral radius of the matrix A^j be less than 1. Recall that the spectral 
radius of a matrix is nothing but its largest eigenvalue. 

III. MULTIVARIATE GENERATING FUNCTION (GF) 

A. Definition for events of first-generation events triggered by a given mother of 
type k 

Among the m types of events, consider the k-th type and its first generation offsprings. 
Let us denote , R^ 2 , ■ ■ ■ , Ri ,m , the number of "daughter" events of first generation of 
type 1,2, ... ,m generated by this "mother" event of type k. With these notations, the 
generating function (GF) of all events of first generation that are triggered by a mother 
event of type k reads 

(5) 

where E [.] represents the statistical average operator. One may rewrite this function in 
probabilistic form 

oo oo m 

A\{ yi ,y 2 , . . . ,y m ) P k( r i,---,r m )Y[y r / , (6) 

n=o r m =0 s=l 
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E 



n 

s=l 



where P fc (ri, . . . , r m ) is the probability that the mother event of type k generates R k > 1 = r\ 
first-generation events of type 1, R k ' 2 = r 2 first-generation events of type 2, and so on. These 
probabilities satisfy to normalizing condition 

oo oo 

^•••^P fe (r 1 ,...,r m ) = l. (7) 

ri=0 r m =0 

The first-order moments or mean values of the numbers of first-generation events of 
different types triggered by a mother of type k are given by 

d 



B. Generating function (GF) for all-generation events triggered by a given mother 
of type k 

The GF A k (yi,y 2 , . . . ,y m ) for all-generation events triggered by a given mother of type 
k is by definition equal to 



A k (y 1 ,y 2 ,...,y m ) := E 



m 

R k,s 



n 



y s 

3=1 



(9) 



where i2 fe,1 ) i? fe ' 2 , . . . , i? fc > m are the numbers of events of all generations and all kinds that are 
triggered by the mother event of type k. 

In order to relate A k (yi, y 2 , . . . , y m ) to A k (yi, y 2 , . . . , y m ), we assume that the first- 
generation daughters can also trigger their own daughters (which are the grand-daughters 
of the initial event) according to the following rules. 

• The numbers of second-generation events that are triggered by each first-generation 
event are statistically independent of the numbers of first-generation events. They 
are also statistically independent of the numbers of second-generation events that are 
triggered by any other first-generation events. 

• Each first-generation event triggers second-generation events according to the same 
laws controlling the triggering of first-generation events by the initial mother event of 
the same type. In other words, the same laws apply to the generation of new events 
from generation to generation, independently of the generation depth. 
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These rules allow us to derive the GF of the numbers of first-generation and of second- 
generation events by performing the following replacement for each variables y q in the ex- 
pression of the GF A\{y x ,y 2 , y m ): 

y q -»■ y q -Al(y 1 ,y 2 ,...,y m ). (10) 

The GF of the numbers of first-generation and of second-generation events that are trig- 
gered by a mother event of type k is given by the following expression in terms of the GF 
A k (yi,y 2 , • • • ,y m ) of the numbers of first-generation events that are triggered by a mother 
event of type k: 

A k (y ll y 2 ,.. . ,y m ) = A\{y t ■ A\(y u y 2 , . . . , y m ), . . . , y m ■ A™(y u y 2 , . . . ,y m )) . (11) 

This equation is valid for all possible values of k = 1, . . . , m. 

By recurrence, one obtain the GF A k +l of the numbers of events of all generations up to 
j + 1 that are triggered by an initial mother event of type k as a function of the GF {Aj} 
of the numbers of events of all generations up to j triggered by an initial mother event of 
type k: 

A* +1 (yi, 2/2, • • • , y m ) = • A](y u y 2 ,...,y m ),...,y m - 4™ (3/1, y 2 , ■ ■ ■ , y m )) ■ (12) 

This equation is valid for all possible values of k — 1, . . . , m and for all possible generation 
levels j = 2 to +00. 

We assume that the above set of recurrence equations ffl2l) for k = 1, . . . , m converges 
to some set of GF's {A k (yi,y 2 , . . . ,y m ); k = l,...,m}. Then, the corresponding GF's 
{A k (yi, y 2 , ... , y m ); k — 1, . . . , m} are solutions of the transcendent equations 

A k (y u y 2 , ...,y m ) = A k (y! ■ A 1 (y u y 2 , . . .,y m ), . . .,y m ■ A m (y 1 ,y 2 , . . .,y m )) , k = 1, . . . ,m . 

(13) 

The equation constitutes the basis for our subsequent analysis. 
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IV. STABILITY CONDITIONS USING MEAN NUMBERS OF EVENTS OF ALL 
GENERATIONS 



A. General relations 

The statistical average of the total numbers of events of type s over all generations that 
are triggered by a mother of type k is given by 

Using ffl3l) . it is straightforward to show that R k,s is solution of 

m 

R h ' s = n ktS + Y,^k,i-R i ' s , (15) 
l=\ 

where n ktS is the mean number of first-generation events of type s triggered by the mother 
event of type k. 

Since expression ( fl5l) holds for all k = 1, . . . , m and s — 1, . . . , m, it can written in matrix 
form 

R = N + NR , (16) 

where N = [ri)~ jS \ is the matrix of the mean numbers of first-generation events and R = [R k,s ] 
is the matrix of the mean numbers of events over all generations. The sum over row indices 
of the elements of the matrix N 

m 

n k = ^2n ktS , (17) 

8=1 

is the mean number of first-generation events of all kinds that are triggered by a mother 
event of type k. 

The solution of the matrix equation (fl6|) is 

N 

i?=-r— -r. (18) 

I — N 

The rest of the paper is concerned with the analysis of particular examples of this general 
solution, worked out for different systems and excitation conditions embodied in different 
forms of the matrix N of the mean numbers of first-generation events. 
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B. Symmetric mutual excitations 



Let us consider the case where 



n k,k — a j 



resulting in the form 



N 



n k , s = b, k ^ s 



abb b 

b a b b 

b b a b 

b b b a 



(19) 



(20) 



for the matrix N of the mean numbers of first-generation events. This form (fT9~|) means that 
events of a given type have identical triggering efficiencies quantified by a to generate first- 
generation events of the same type. They also have identical efficiencies quantified by b to 
trigger first-generation events of a different type. In other words, the mean number of first- 
generation events of type k triggered by a mother event of the same type k is independent 
of k. And the mean number of first-generation events of any type s 7^ k triggered by any 
another event of a different type k is independent of k and s. 

As a consequence, the mean number of first-generation events of all kinds that are trig- 
gered by a mother event of some type k, as given by ffTTl) . is independent of k and given 
by 

rtj, = n = a + (m — 1)6 
It is convenient to introduce the factor 

b 

Q = ~ 

a 

comparing the inter-types with the intra-type triggering efficiencies. Using definition ( 122]) 
and equality ( l2Tj) . we obtain 

n nq 



for all k . 



(21) 



(22) 



1 + (m — l)q ' 1 + (m — l)q 

Two limiting cases are worth mentioning: q = (independent types) and q 
equivalent types): 



n, 



b\ q=0 = 0, 



1 9=1 



b\ 



q=l 



n 
rn 



(23) 
1 (fully 

(24) 
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The solution ( ITS]) implies that the matrix R possesses the same structure as the matrix 
N, with identical diagonal elements R k,k and identical off-diagonal elements R k,s (for k 7^ s), 
given respectively by 

£>fc,fc = ^ ) k = l,...,m, 

1 — n l + n(q — l) + q(m—l) , nr , 

(25) 

f)k,s _ n <7 7 / 

XL — " J ~ ~ ~ . rC t^- 5. 

1 — n 1 + n{q — 1) + g(m — 1) 

Therefore, the mean R k of the total number of events of all kinds, that are triggered by 
some given mother event of a given type k, is given by 

m 

R k = J2R k ' s = -=R, Vg . (26) 

s=l 

This expression has a simple interpretation, resulting from equation ( 12~T|) and the process 
of triggering. Indeed, by definition of n in (T2T|) . there are on average n first-generation 
events of all kinds that are triggered by a mother event of some type k. Each of these 
first-generation event triggers on average n second-generation events of all kinds, leading 
to a total contribution n 2 for the number of second-generation events. Counting all the 
generation cascades, we obtain n + n 2 + n 3 + which is nothing but the result ( 1261) . 

As for the mono-variate Hawkes process, the dynamics is stable (sub-critical) for n < 1 
and unstable (super-critical or exponentially explosive) for n > 1. As usual, the critical point 
occurs when there is exactly n = 1 first-generation events of all kinds that are triggered 
by a mother event of any type. There is not qualitative difference between this multi- 
variate Hawkes process with the structure (f20i) of mutual excitations and a mono-variate 
Hawkes process, once the branching ratio n defined as the average number of first-generation 
daughters from a given mother is generalized into its natural extension fT2TT) . 

C. Two-dimensional mutually and self-excited Hawkes process 

With only two types of events, a detailed analysis can be performed, with the discovery 
of new qualitative regimes. 
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1. Stability analysis 



With two types of events, the 2x2 matrix N of the mean numbers of of first-generation 
events can be kept fully general and is noted as 



N 



niA ni )2 
n 2 ,i n 2 , 2 



(27) 



(28) 



The solution (TT8]) is a 2 x 2 matrix R with elements R ,s given by 

- ll m t i + 711,2712,1 - "1,1^2,2 519 %,2 

-M. — — , K ' — — , 

V V 

d2 1 ^2,1 ^2 2 %,2 + 711,2712,1 - 711,1712,2 
/i — — , it ' — — , 

V V 

where 

V = 1 + 711,17^2,2 - 711,2712,1 - 77-1,1 - 772,2- (29) 

In order to determine the stability of two-dimensional mutually and self-excited Hawkes 
process, we study the mean numbers n\ and n 2 of first-generation events triggered by a 
mother event of the first and second kind, respectively. They are given by the sums of the 
two row elements of the matrix N: 

771 = 771,1 + 7-71,2, " 2 = " 2 ,1 + 77 2 , 2 . (30) 

It is convenient to use a representation of the elements of the matrix N similar to (1231) : 

ni 77igi 77 2 g 2 77 2 f . 

n l,l = T~ > n l,2 = ~r, > "2,1 = , 77 2 , 2 = — . (31) 

l + 1 + gi 1 + 92 1 + 92 

As in (I22p . 91 and 92 quantify the relative strengths of inter-type compared with the intra- 
type triggering efficiencies: 91 = 77,1, 2/77.1,1 and 92 = "2,1/^2,2- The limit 91 = q 2 = reduces 
to two independent self-excited Hawkes processes. 

The solutions R k,s given by ( 128]) are finite as long as the spectral radius A(ni,n 2 ) of 
the matrix iV remains smaller than 1. This defines the sub-critical regime. The system 
becomes critical (respectively super-critical) when the spectral radius A(ni,n 2 ) is equal to 1 
(respectively larger than 1). One can show that the set {n\\ n 2 ) such that the denominator 
V given by (f2"9"|) is identically zero is critical, i.e., corresponds to a unit spectral radius, if V 
remains positive in the domain bounded by the semi-axes [771 G (0; +oo), 77 2 £ (0; +00)] and 
the curve T>{rii] n 2 ) = 0. 
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In order to study the three regimes (sub-critical, crtical and super-critical), it is convenient 
to express V(ni,n 2 ) as a function of n\ and n 2 , using (1311 : 




Fig. 1: Critical line T>(ni,n 2 ) = and three domains in the plane (ni,n 2 ): 1) 
subcritical, where both n\ and n 2 are smaller than 1; 2) quasi subcritical, where one 
of the mean numbers (n\, n 2 ) of first-generation events is larger than one but the mean 
numbers of events of all generations are finite; 3) supercritical region, where all mean 
numbers R k,s are infinite. 

The critical line V(ni,n 2 ) = is shown in figure 1 in the plane (ni,n 2 ), together with 
three domains. 

1. For ni < 1 and n 2 < 1, V(rii,n 2 ) > and the system is subcritical. The conditions 
ri\ < 1 and n 2 < 1 mean that the cascade of events over all generations do not blow 
up for each of the two types of event triggering. 

2. The domain indicated in figure 1 as "quasi subcritical" is such that one of the mean 
numbers (ni,n 2 ) of first-generation events is larger than one but the mean numbers 
of events over all generations remain finite since T>(rii, n 2 ) > 0. Intuitively, the super- 
critical regime of one of the event types is damped out by the triggering of the second 
type of events which is subcritical. The two extreme boundaries {rii — 1 + qi, n 2 — 0) 
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and (ni = 0; n 2 = 1 + q 2 ) exemplify this point as they correspond respectively to 
72i ; i = 1; ri2 = and n\ = 0; 712,2 = 1. 

3. In the domain indicated at "supercritical" in figure 1, the mean number of events of 
both types summed over all generations goes to infinity. This occurs of course if both 
ni and n 2 are larger than 1 but also when one of them is smaller than 1 if the other 
one is sufficiently large. In this later case, the damping offered by the second type is 
not sufficient to stabilize the triggering process. This is the runaway explosive regime. 

4. The downward sloping line defines the critical domain T>(n 1 ,n 2 ) = separating the 
quasi subcritical and the supercritical regimes. 

As the critical line 'D{n\^n 2 ) = is approached from within the subcritical regime, the 
total mean number R 1 = R 1,1 + R 1 ' of events of all types and over all generations, 

ni(l + qi + g 2 + <?i<?2 - n 2 (l - gift)) 



R^R^+R 1 ' 1 



(33) 



^1 + g 2 )(l + qi~ ni) - (1 + q\ - n x {\ - qiq 2 ))n 2 ' 

grows to finally diverge on the line, as shown in figure 2 for a particular example. 

In this example, q% = 0.2, q 2 = 0.4, ni = 0.8. The value qi/(l +qi) = 16.7% is the fraction 
of first-generation events generated by a mother of the first type which are of the second 
type. The value q 2 /{l + q 2 ) = 28.6% is the fraction of first-generation events generated by 
a mother of the second type which are of the first type. 
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Fig. 2: Linear- log plot of the mean number R 1 of events of all types and over all 
generations given by expression (|33p as a function of ri2 for q\ = 0.2, c/2 = 0.4 and 
m = 0.8. Note that R 1 remains finite even when becomes larger than 1 up to 
a critical value n\ = 1.206896... for which T>(ni,ri2) = at which it diverges. For 
ri2 = 1.2 for instance, R 1 = 144. 

2. Strong asymmetry in mutual triggering 

It is instructive to consider the limiting case where one type of events triggers many more 
events of the other type than the reverse. Mathematically, this corresponds to 



which means that the fraction of first-generation events generated by a mother of the first 
type which are of the second type is much larger than the fraction of first-generation events 
generated by a mother of the second type which are of the first type. In this case, the critical 
line T>(ni,n2) = becomes almost rectangular, as illustrated in figure 3. 



(34) 



1.2r 




0.4 



0.2 - 
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1.5 
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Fig. 3: Plots of the critical line T>{n\,n2) = for the four following values of (q\; q2): 



(0.1;0.01), (0.5;0.01), (1;0.01), (1.5;0.01) from left to right. 
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Let us consider in more details the limiting case q\ = q while q2 = 0, for which the 
relations fl3TT) transform into 



ni.i = — — , ni j2 = — — , n 2 ,i = , n 2j2 = n 2 . (35) 
l + q 1 + q 



As a result, the relations (I28I) and (|29|) become 



l + q- (l + ?-ni)(l-n 2 )' ^ 

i? 2 - 2 = i? 2 ' 1 = 0. 

1 — n 2 

The finiteness (subcritical behavior) of the number Rr' 1 is controlled solely by rii, which 
must be smaller than l + q. As n\ tends to 1 + q, Br' 1 goes to infinity, expressing the 
transition to the supercritical regime. In contrast, there are two mechanisms leading to the 
divergence of B 1,2 . 

1. As rii tends to 1 + q, R 1 ' 2 goes to infinity, as the number of events of the first type 
itself diverges, each of these events triggering a significant fraction of events of the 
second type, at each generation. In other words, the divergence of R 1 ' 2 is controlled 
by or slaved to that of R 1 ' 1 , which reflects the triggering efficiency of events of type 
one. 

2. The number B 1,2 of events of the second type generated over all generations by a 
mother event of type one diverges when n 2 — > 1, even if ni < 1 + q. A mother 
event of type one triggers events of type two at each generation, each of these events 
only triggering events of their own kind. Thus, the number of events of the second 
type diverges when the self-triggering parameter (or "branching ratio" ) n 2 reaches its 
critical value n\ = 1. 

This implies in particular that, when q or n 2 are sufficiently large, the following inequality 
holds: B 1,2 > R 1 ' 1 . The general condition for this to be true is n 2 + q > 1. Intuitively, for a 
fixed self-triggering ability n 2 , there must be sufficiently many events of type two generated 
by events of type one: q > 1 — ri 2 . Alternatively, for a fixed fraction q/(l + q) of first- 
generation events of type two generated by events of type one, the branching ratio of type 
two events must be sufficient large: n 2 > 1 — q. 

Note that, if both n 12 and n 2 \ are positive, corresponding to a non-degenerate case, then 
events of any type do trigger events of the other type. As a consequence, either R 1 and R 2 
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are both finite or both infinite. This is not the case for independent or semi-independent 
systems. A system is independent if n\ 2 = 7121 = 0, i.e., events of different types live their 
separate "lives" without any inter-mutual triggering. In this case, for instance, if n\ < 1, 
while n% > 1, then the events of the first type form a subcritical set, while the events of the 
second type form a supercritical system. Two systems are semi-independent if only one of 
the two cross-terms n\ 2 and 7121 is equal to zero. Suppose for instance that n 2 \ = 0. Then, 
if Ft}' 1 is finite, then R 1 ' 2 might be finite or infinite, just because events of the second type 
cannot trigger events of the second type and an infinite value of R 1,2 remains compatible 
with a finite value of R 1 ' 1 . In contrast, if R 1 ' 1 is infinite, then R 1 ' 2 is necessarily infinite, 
because the infinite number of events of the first type trigger an infinite number of events of 
the second type, since n i2 7^ 0, even if self-triggering of events of type two is zero {n 22 = 0). 
Moreover, a main peculiarity of degenerate (independent or semi- independent) systems is a 
strongly rectangular critical curve. 




0.2 0.4 0.6 0.8 1 1.2 

ni 



Fig. 4: (color online) Linear-log dependence of the mean numbers R ' and R > 2 given 
by (f36j) as a function of m, for q = 0.2 and for two values n 2 = 0.5 and n 2 = 0.9 

Figure 4 shows in linear-log scale the number R 1,1 (respectively R 1 ' 2 ) of events of the first 
type (respectively second type) over all generations that are triggered by a mother event of 
the first type, given by f l5B|) . as functions of ni, for q = 0.2 and for two values n 2 = 0.5 and 
n 2 = 0.9. One can observe that, for n 2 = 0.9, R 1,2 is larger than R 1,1 for any n\ 6 (0, 1.2). 
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D. One-dimensional chain of directed triggering 



The case of a strong asymmetry in mutual triggering discussed in subsection IIVC2I for 
the two-dimensional case can be generalized to the case of m > 2 different event types. We 
consider a chain of directed influences k — > k + 1 where the events of type k trigger events 
of both types k and k + 1 only, and this for k = 1, 2, m. This is captured by a form of 
the matrix N which has only the diagonal and the line above the diagonal with non-zero 
elements. 

As the simplest example, consider first the matrix N 



N 



X s 0... 

x s 0... 

x s 0... 

x s 0... 



(37) 



where 



n 



nq 



X = t— , 8 = -f- . (38) 
1 + q 1 + q 

Reduced to a two-dimensional system m — 2, this corresponds to the particular case of 
subsection II VC 21 for which 77-2,2 = ^2 = + q) in the notations of expressions fl35|) . 

because the diagonal elements are taken all equal. For m = 3, a financial example is that 
the events of the first type correspond to fundamental news, the events of the second type 
are the price jumps of a leading market such as the US (assuming no feedbacks of prices 
on news) and the events of the third type correspond to the price jumps of a secondary 
market, such as the Russian stock market (assuming to effect of the Russian market on the 
US market). 

Assuming that the mother event is of type k, the solution of equation ffl8|) for this case 
(|37j) with flU is given by R s > s = R k ' s = R s > k = for 1 < s < k and 

X n 



R 



k.k 



^s—k 



R 



k,s 



i-x 

(1+?)- 



1 + q — n 

(nq) s ~ k 



(39) 



\s-k+l ' 



S > k ■ 



The fact that the critical point is solely controlled by a single critical value n c — 1/(1 + q) 
results from the fact that all diagonal elements of the matrix fl37|) are equal. 
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E. One-dimensional chain of nearest-neighbor triggering 

A natural extension to the above one-dimensional chain of directed triggering discussed 
in section IIVDI is to include the possibility of feedbacks from events of type k + 1 to type 
k. The simple example is to consider symmetry mutual excitations confined to nearest 
neighbors in the sense of event types: k k + 1. Mathematically, this is described by a 
symmetric matrix N of the average numbers jS of first-generation events of different types 
triggered by a mother of a fixed type. Figure 5 provides the geometrical sense of matrix iV 
(HO]) for m — 6, where the circles represent the six types of events and the arrows denote 
their mutual excitation influences. 




Fig. 5: Geometric sense of the matrix ./V for a one-dimensional chain of nearest- 
neighbor triggering. 

Here, we consider that all diagonal elements are equal to some constant \ (same self- 
triggering abilities) and all off-diagnoal elements are equal to some different constant s 
(same mutual triggering abilities). The elements ni jm and n mj i are also equal to s to close 
the chain of mutual excitations between events of type 1 and of type m. Restricting to 
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m = 6 for illustration purpose, the corresponding matrix N reads 

X s s 

s x s 

s x s 

s x s 

s x s 

s s x 



N 



(40) 



where 



n 



X 



nq 



l + q 



X + 2s 



n 



(41) 



2(1 + 9) 

As before, the parameter q quantifies the "strength" of the interactions between events of 
different types. Here, n represents the total number of first-generation events of all types 
that are generated by a given mother of fixed arbitrary type. 

The solution of equation ( fl8l) for this case is given by a circulant structure 

A B C D C B 

B A B C D C 

C B A B C D 

D C B A B C 

C D C B A B 

B C D C B A 



R 



(42) 



where 



4n(l — n + q) 3 — (1 — n + q)(5n — 2q — 2)n 2 q 2 — n 4 g 4 
4(1 - n + qY - 5(1 - n + q) 2 n 2 q 2 + n^q* 

nq(l + g)(2(l — n + q) 2 — n 2 q 2 ) 



A(n,q) 



(43) 



B(n,q) 
C(n } q) 
D(n,q) 

It is straightforward to check that the common denominator 4(1 — n + g) 4 — 5(1 — n + 
q) 2 n 2 q 2 + n A q A to these four numbers A(n, q), B(n, q), C(n, q) and D(n, q) 



4(1 


- n + q) A 


-5(l-n 


_l_ qYn 2 q 2 


+ n 4 q 4 




n 2 q 2 


(1 + <?)(! - 


-n + q) 




4(1 


-n + qY 


-5(1 -n 


_l_ qYn 2 q 2 


+ n 4 q 4 






n 3 q 3 (l + 


q) 




4(1 


-n + qY 


-5(l-n 


_l_ qYn 2 q 2 


+ n 4 q* 



1. does not vanish for any q values for n < 1, 



20 



2. vanishes for any q values at n = 1, and 

3. can vanish at up to four values of q for n > 1. 

This means that the system resulting from the structure of mutual excitations between 
different event types represented by the matrix (142!) is always in the subcritical regime for 
n < 1 and becomes critical at n = 1 as for the decoupled or monovariate case. 
The behavior of the ratios 



U2 -R^-A> **-WA' n '-R^~A 



(44) 



is shown in figure 6 as a function of n for q = 0.1. Here, the mother event is of type 1 and 
the curves illustrate the progressing dampening of the cascade of triggering proceeding from 
type to type via nearest-neighbor mutual excitations. 




0.94 0.96 

n 



Fig. 6: (color online) Top to bottom: ratios IZ2, ^3 and IZ4 defined in (|44p of the 
mean numbers of events of all generations generated by a mother event of type 1 for 
q = 0.1. 



F. Subcriticality measure in the case of two types of events 

For a system with two types of events, the explicit relation defining the critical curve is 

n 2 = G(ni,q 1 ,q 2 ) = j— K r f. 7 , ^ m < 1 + q u (45) 

(1 + gi) - (1 - qmjni 
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which is represented in figure 7. One can verify that the point (ni = l,ri 2 = 1) is always on 
the critical curves. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

m 



Fig. 7: (color online) Plots of critical curves n2(ni) for systems with two types 
of events, for the following pairs of parameters: bottom up on the left side of the 
curves, we have {q\ = 0,(72 = 0), (gi = 0A,q 2 = 0.6), {q\ = 0.7,^2 = 0.2) and for 
(gi = 0.8,g 2 = 1.1). 

One can observe that, for qi > and q 2 > 0, the subcritical domain is significantly larger 
than for systems in which event types are independent, i.e., do not mutually trigger each 
other, corresponding to qi = q 2 = 0. It is illuminating to introduce a quantitative measure 
of the domain of subcriticality, here chosen for the two-dimensional case as the surface S of 
the subcritical domain: 

rl+q 

S{qi,q2)= / G(n 1 ,qi,q 2 )dn 1 . (46) 
Jo 

For independent types (gi = q 2 = 0), 5(0, 0) = 1. The general expression of the subcriticality 
measure S(qi,q 2 ) is obtained as 

cr \ ft , \ft 1 x 1 - 9i?2 [1 - ln(gig 2 )] 

S{qu<h) = (l + gi)(l + «0 t, ^ • 47 

(1 -qiq 2 ) 2 

In the particular case of a chain of directed triggering in the space of event types studied in 
subsection IIVDI for which the critical curve is rectangular, we find 

S(q) := lim S(q,q 2 ) = l + q. (48) 
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In the case of symmetric triggering among different event types (qi = q2 = q), we obtain 

1 — q 2 + q 2 ln(g 2 ) 



S(q,q) 



(49) 



(1-q) 2 

These last two functions S(q) and S(q, q) are depicted in figure 8. The main insight is that 
mutual triggering tends to provide a significant extension of the stability domain. 




Fig. 8: (color online) Dependence of the subcriticality measures S{q) and S(q,q) 
given respectively by expression (148 p (straight line) and (|49|) (concave curve). 



V. CONCLUSION 



Considering the class of multivariate self-excited Hawkes point processes, we have pre- 
sented the general theory of multivariate generating functions to derive the number of events 
over all generations of various types that are triggered by a mother event of a given type. 
This has allowed us to discuss in details the stability domains of various systems, as a func- 
tion of the topological structure of the mutual excitations across different event types. In 
particular, we have studied the case of symmetric mutual excitation abilities between events 
of different types, the case of just two different types of events, the case of a one-dimensional 
chain of directed triggering in the space of event types and the case of a one- dimensional 
chain in the space of event types with nearest-neighbor triggering. The main insight is that 
mutual triggering tends to provide a significant extension of the stability (or subcritical) 
domain compared with the case where event types are decoupled, that is, when an event of 
a given type can only trigger events of the type. 
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